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Abstract. Far-infrared spectra (43-197 /im) of 34 nearby galaxies obtained by the Long Wavelength Spectrometer 
(LWS) aboard the Infrared Space Observatory (ISO) were analyzed to investigate the general properties of in- 
terstellar matter in galaxies. The present sample includes not only normal galaxies but also starbursts and 
active galactic nuclei (AGNs). Far-infrared forbidden lines, such as [Cn]158/im, [Ol]63/im, [Nil] 122 /im, and 
[Oni]88/^m, were detected in most of the sample galaxies. [Ol]145^m line was detected in 13 galaxies. The line 
fluxes of [Cll]158/xm and [Nil] 122 /im relative to the total far-infrared flux (FIR) decrease as the far- infrared 
color becomes bluer, while the ratio of the [Ol]63/im flux to FIR does not show a systematic trend with the 
color. The [Olll]88/im to FIR ratio shows a large scatter with a weak trend of increase with the color. AGNs 
do not show any distinguishable trend from normal and starburst galaxies in the far-infrared spectra, suggesting 
that the far-infrared emission is mainly driven by star-formation activities even in AGNs. We estimate the phys- 
ical conditions of photodissociation regions (PDRs) in the sample galaxies, such as the far-ultraviolet radiation 
field intensity Go and the gas density n by assuming that all the observed [Ol]63/^m and far- infrared contin- 
uum emissions come from PDRs. Comparison with PDR models indicates that Go ranges from 10^ — 10^ and 
n~10^-10*c] 
to CO (J = 1 

[Cll]158/im emission originates from PDRs and attribute the rest to the emission as coming from low-density 
diffuse ionized gas. The estimated intensity of [Cll]158/im from the ionized gas is compatible with the observed 
intensity of [N ll] 122 fim if both lines come from the same diffuse ionized gas. The present analysis suggests that 
the decrease in [C ll] 158 fim/FIR with the far-infrared color may not be accounted for by the decrease in the 
photoelectric heating efficiency owing to the increase in positive charges of dust grains because a measure of the 
efficiency, Go/n, is found to stay constant with the far-infrared color. Instead the decrease can be interpreted in 
terms of either the increase in the coUisional de-excitation of the [C ll] transition due to the increase in the gas 
density or the decrease in the ionized component relative to the far-infrared intensity suggested by the decrease 
in [N ll] 122 fim/FIR. Based on the present analysis, we derive average relations of the far- infrared color with Go 
and n in galaxies, which can be applied to the investigation of interstellar matter in distant galaxies. 

Key words, galaxies: ISM - infrared: ISM - lines and bands: ISM - radiation mechanisms: thermal 



'^ . The present results also suggest that n varies proportionally with Go. The ratio of [C ll] 158/im 
0) line emission supports the linear increase in n with Gq. We estimate that about a half of 



1. Introduction 

Mechanisms for heating interstellar gas include colli- 
sions, radiation from stars, shocks, and cosmic rays. 
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Examination of the spectral lines that cool the gas can 
help determine the dominant excitation mechanisms and 
conditions. Studies of particular regions in our Galaxy and 
observations of external galaxies have suggested that stel- 
lar ultraviolet radiation can ionize vast volumes of a galaxy 
and that the far-ultraviolet (FUV) radiation impinging on 
neutral cloud surfaces is responsible for a large fraction of 
the observed far-infrared (FIR) spectral line emission that 
cools the gas (Cravifford et al. 1985; Stacey et al 
l99l| ) 



Shibai et al. 



Tielens & HoUenbach ( [1985D defined 
photodissociation regions (PDRs) as "a neutral region 



1991 
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where FUV radiation dominates the heating and/or some 
important aspect of the chemistry" . Thus PDRs include 
most of the atomic gas in a galaxy, both in diffuse clouds 
and in dense regions (for a recent review, HoUenbach & 
Tielens p^97| ). 

[Cii]158/ini and [Oi]63/^ni lines are important 
coolants in PDRs, while gas heating in PDRs is thought 
to be dominated by energetic photoelectrons ejected from 
dust grains fo llowin g FUV photon absorption (Watson 



1972| ; de Jong |1977|) . For galactic nuclei and star- forming 



regions in the spiral arms, most of the observed [C ii] line 
emission arises from PDRs on molecular cloud surfaces. 
However, integrated over the disks of spiral galaxies, a 
substantial fraction may also arise from "standard" atomic 
clouds, i.e., th e cold neutral medi um ga s regions (CNM) 
(Madden et al. 1993 ; Bennett et al. 1994 ) or from extended 
low-d ensity warm ionized gas regions (ELDWIM) (Heiles 
1994). Contributions from the various gas phases can be 



estimated by obser vations of seve ral FI R forbidden lines 
(Carral et al. 1994 Luhman et al. 



1998) 



With the LWS on board the ISO (Clegg et al. |1996D 
it now becomes possible to measure far-infrared lines 
not only from infrared-bright galaxies but also from nor- 
mal galaxies (for latest reviews, Genzel & Cesarsky 2000; 
Fischer 2000). Lord et al. (1996) made observations of sev- 
eral FIR fine-structure forbidden lines in a normal galaxy 
NGC5713 and found a fairly high [ Cu]15 8^lUi to the FIR 
intensity ratio. Smith & Madden (1997) made observa- 
tions of [Nii]122/^m and [Cii]158/im lines for five spi- 
ral galaxies in the Virgo cluster and found enhanced ra- 
tios of [Cii] to CO(J =1 — 0) intensities in two of the 
five galaxies. They interpreted the results in terms of ei- 
ther low metallicity or an increa se in the contribution 
from the CNM. Braine & Hughes (1999) obtained a com- 
plete FIR spectrum of a normal disk galaxy NGC4414 
and investigated the physi cal co nditions of the interstel- 
lar medium. Leech et al. ( 1999| ) presented observations 
of [Cii] 158 /im in 19 quiescent galaxies in the Virgo clus- 
ter. They indicated a good correlation between [C ii] and 
far-infrared intensities and a trend of increasing [C ii]-to- 
far-infrared intensity ratio with increasing galaxy lateness, 
which has been shown to be related to th e star- formation 
rate (Pierini et al. 1999| ). Malhotra et al. ( 2001a ) reported 
observations of 4 early-type galaxies with ISO/LWS and 
interpreted the observed low ratio of [Cii]158/im to far- 
infrared intensities in terms of the soft radiation field in 
the target galaxies. Hunter et al. ( 200l| ) presented obser- 
vations of 5 irregular galaxies with ISO/LWS and found 
strong [Cii] 158 /xm emission relative to the far-infrared 
continuum. 

Malhotra et al. ( |1997| , |20011tD have investigated the 
far-infrared properties of 60 nearby normal galaxies based 
on line-spectroscopic observations mainly of [Cii]158/im 
and [Oi]63/im line emissions. They complemented their 
line data with IRAS photometry to estimate the far- 



Several interpretations have been proposed for the trend 
(M alhot ra et al. |1997| ; Ge nzel fc Cesarsky |2000| ; Helou et 
al. 2001). Malhotra et al. (2001b) favor the interpretation 
of the decrease in the photoelectron yield owing to the 
increase in positive charges of dust grains under strong 
ultraviolet radiation. 

The FIR continuum emission shorter than 60 fiia is 
dominated by the emission from stochastically heated very 
small grains (Des ert et al. 1986 ; Dweket al. 1997 Onaka 
^OOOt Dale et al. ^00 1| ; Li & Draine |200l| ). The spectral 
energy distribution (SED) longer than 100 /im is crucial 
to correctly estimate the thermal emission from submi- 
cron dust grains and understand the FIR SED of galax- 
ies. In this paper we investigate the FIR properties of 
nearby galaxies based on LWS full grating spectra from 
43 to 197 /im including both line and continuum emission. 
The continuum spectra longer than 100 /xm enable better 
estimates of the average temperature of submicron dust 
grains as well as the strength of the interstellar radiation 
field. Because the aperture size of LWS is large {^ 80 "), 
the spectra of galaxies include contributions from various 
interstellar regions within the galaxies. This paper inves- 
tigates the mechanism of gas heating with the aim of a 
better understanding of the global physical conditions of 
the interstellar medium in galaxies. 



2. Data reduction and observational results 

Full grating scan spectra of 43 — 197 /im were obtained 
for 9 galaxies with the LWSOl mode in the open time 
programs of TONAKA.IRGAL and GALIR. In addition, 
we extracted LWSOl full grating scan data of 25 galaxies 
from the ISO archival database for a total of 34 nearby 
galaxies LWS spectra that were analyzed in the present 
study. The sample includes various types of galaxies, rang- 
ing from active galactic nuclei (AGNs), starburst, to nor- 
mal galaxies. Table |l| lists the present sample, where the 
flux density at 60 /<m and the FIR color i?(60/100) are 
derived from the LWS spectra convolved with the IRAS 
band filters (for the data reduction, see below). We ex- 
cluded Arp220 and II ZW40 from the present analysis, 
both of which have also been observed with LWSOl. The 
former shows a spectrum optically thick even in the far- 



infrared continuum intensity FIR (Helou et al. 1988) and 
found a trend that the ratio of the [C ii] line intensity to 
FIR decreases with the far-infrared color becoming bluer. 



infrared (Fischer |2000[ Malhotra et al. |20011| ), for which 
the present simple analysis cannot be applied, and the 
continuum of the latter object is faint and suffered from 
the uncertainty in the dark current estimate. The spectral 
resolution (A/AA) of LWOl mode was about 200 and the 
beam size of LWS was estimated to be 60" — 80" (Cry 
EOOO). The sample galaxies have optical sizes of 1' — 30' 
and thus for some galaxies only the central portion was in- 
cluded in the LWS beam. Smith & Harvey ( 1996| ) reported 
observations of far-infrared emitting regions in external 
galaxies from the Kuiper Airborne Observatory (KAO), 
indicating that most far-infrared emission is concentrated 
in the central 30" regions of the galaxies. LWS observa- 
tions are supposed to detect most of the far-infrared emis- 
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. List of the galaxies in the present analysis 








Galaxy 


Type 


Morphology 


^0(60)" [Jy] 


R(60/100)'' 


visual size '] 


Gen A 


AGN 


SO pec, Sy2 


98.4 


0.43 


25.7 X 20.0 


Gircinus 


AGN 


SA(s)b:, Sy2 


334.1 


0.71 


6.9 X 3.0 


IC2554 


normal 


SB(s)bc pec: 


15.9 


0.52 


3.1 X 1.3 


IRAS00506+7248 


normal 




24.1 


0.66 




IRAS13242-5713 


normal 


S... 


89.2 


0.70 


1.1 : xO.2 


M51 


normal 


SA(s)bc pec, HIT, Sy2.5 


36.4 


0.38 


11.2 x 6.9 


M82 


starburst 


10, Sbrst, HII 


1486.6 


0.87 


11.2 X 4.3 


M83 


starburst 


SAB(s)c, HII, Sbrst 


138.0 


0.60 


12.9 X 11.5 


Maffei2 


normal 


SAB(rs)bc: 


94.8 


0.47 


5.82 X 1.57 


NGG1068 


AGN 


(R)SA(rs)b, Syl, Sy2 


206.3 


0.62 


7.1 X 6.0 


NGG1097 


starburst 


(R'i:)SB(r'l)b, Syl 


49.6 


0.54 


9.3 X 6.3 


NGG1365 


AGN 


(R')SBb(s)b, Syl.8 


92.0 


0.52 


11.2 X 6.2 


NGG2146 


starburst 


SB(s)ab pec, HII 


163.3 


0.65 


6.0 X 3.4 


NGG253 


starburst 


SAB(s)c, HII, Sbrst 


1044.7 


0.70 


27.5 X 6.8 


NGG3256 


starburst 


Pec, merger, HII 


107.8 


0.71 


3.8 X 2.1 


NGG3690 


starburst 


IBm pec, HII 


121.7 


1.00 


2.0 X 1.5 


NGG4038 


starburst 


SB(s)m pec 


21.5 


0.52 


5.2 X 3.1 


NGG4041 


normal 


SA(rs)bc: 


11.7 


0.43 


2.7 X 2.5 


NGG4414 


normal 


SA(rs)c? 


24.8 


0.38 


3.6 X 2.0 


NGG4945 


starburst 


SB(s)cd: sp, Sy2 


577.2 


0.50 


20.0 X 3.8 


NGG520 


starburst 




37.9 


0.65 


1.9 xO.7 


NGG5430 


starburst 


SB(s)b, HII, Sbrst 


9.2 


0.51 


2.2 X 1.1 


NGG5937 


normal 


(R')SAB(rs)b pec 


10.0 


0.46 


1.9 X 1.1 


NGG6156 


normal 


(R'i)SAB(rs)c 


20.9 


0.61 


1.6 X 1.4 


NGG6240 


starburst 


10: pec, LINER, Sy2 


25.8 


0.87 


2.1 X 1.1 


NGG6764 


starburst 


SB(s)bc, LINER, Sy2 


4.6 


0.50 


2.3 X 1.3 


NGG6810 


normal 


SA(s)ab:sp, Sy2 


17.0 


0.48 


3.2 X 0.9 


NGG6824 


normal 


SA(s)b: 


6.8 


0.43 


1.7 X 1.2 


NGG6946 


starburst 


SAB(rs)cd, HII 


58.6 


0.50 


11.5 X 9.8 


NGG7469 


starburst 


(R')SAB(rs)a, Syl.2 


24.7 


0.69 


1.5 X 1.1 


NGG7552 


starburst 


(R')SB(s)ab, HII, LINER 


74.7 


0.59 


3.4 X 2.7 


NGG7582 


AGN 


(R'i)SB(s)ab, Sy2 


51.9 


0.64 


5.0 X 2.1 


NGG7714 


starburst 


SB(s)b:pec, HH, LINER 


12.7 


0.97 


1.9 X 1.4 


NGG891 


normal 


SA(s)b? sp, HH 


24.3 


0.31 


13.5 X 2.5 



a Flux density at 60 /^m derived from the LWS spectra convolved with the IRAS band filters. 

b Ratio of the 60 /xm and 100 /xm flux densities (in Jy) derived from the LWS spectra with the IRAS band filters. 



sion and therefore probe the properties of central part of 
the galaxies. 

In the present study, we used the Standard Processed 
Data (SPD) of off-line processing (OLP) version 9 prod- 
ucts provided by the ISO data center. The dark current 
and the drift in the detector responsivity were corrected by 
using the LWS Interactive Analysis software (LIA version 
7.3)ffi The ISO Spectral Analysis Package (ISAP version 
1.6a)n was then used for further data reduction. The con- 
tinuum spectra were stitched together by shifting each de- 
tector signal with the offset method, adjusted to the SW5 
channel in most galaxies. In some cases where the SW5 
channel is noisy, the adjustment was made referring to 



^ LIA is a joint development of the ISO-LWS instrument 
team at RAL, UK, the PI institute and IPAC. 

^ The ISO Spectral Analysis Package (ISAP) is a joint de- 
velopment by the LWS and SWS Instrument Teams and Data 
Centers. Gontributing institutes are CESR, IAS, IPAC, MPE, 
RAL and SRON. 



the LWS or LW4 channels. The offsets between the detec- 
tors were typically less than 20%. The line flux, the total 
far-infrared flux, and the dust temperature of the contin- 
uum emission were derived by ISAP. The 60 and 100 /im 
flux densities from the LWS spectra were found to agree 
with the IRAS data within about 20%. The uncertainties 
in LWS spectra were suggested to be about 15 — 20% in 
previous works (e.g., Braine and Hughes 1999| ; Unger et 
al. 2000) and we adopt 20% errors for the flux uncertainty. 



The continuum emission shorter than the 60 /zm region 
is affected by the contribution from very small grains. To 
derive a typical temperature of submicron grains Td in 
each galaxy, we fitted the LWS spectrum for A > 80 /im 
with the following equation: 



F(A) = O To. 



55 



0.55/im\ 



BxiTi), 



(1) 



Table 2. Summary of the observational results 



Name 


FIR'' 

(lO-inVm-2) 


(K) 


[0...](52)<= 
(10-"Wm-2) 


[Oi](63)<' 
(10-"Wm-2) 


[0 III] (88)'= 
(lO-i^Wm-2) 


[Nii](122) ^ 
(10-"Wm-2) 


[Oi(145)]» 
(lO-i^Wm-2) 


[Cii](158)'' 
(10-"Wm-2) 




Cen A 


9.1 ± 2.0 


35.2 ± 1.5 


4.4 ± 1.2 


19.2 ± 3.9 


6.3 ± 1.4 




1.0 ± 0.2 


29.0 ± 5.8 




Circinus 


22.3 ± 3.0 


42.4 ± 1.2 


8.6 ± 2.6 


23.0 ± 4.6 


8.5 ± 1.8 




1.8 ± 0.5 


26.5 ± 5.3 




1C2554 


1.2 ± 0.4 


35.1 ± 2.3 




2.3 ± 0.5 


1.3 ± 0.3 






3.6 ± 0.7 




1RAS00506+7248 


1.5 ± 0.5 


39.1 ± 2.8 




2.7 ± 0.5 


1.9 ± 0.4 






2.9 ± 0.6 




IRAS13242-5713 


5.8 ± 2.2 


40.9 ± 3.1 




4.0 ± 0.8 


3.2 ± 0.7 






5.3 ± 1.1 




M51 


3.5 ± 0.6 


29.9 ± 1.1 




4.4 ± 0.9 


0.8 ± 0.3 


2.1 ± 0.4 




10.4 ± 2.1 




M82 


101.2 ± 22.3 


49.8 ± 2.2 


102.5 ± 21.1 


169.4 ± 33.8 


90.6 ± 18.2 


20.5 ± 4.2 


14.6 ± 2.9 


128.0 ± 25.6 




M83 


9.6 ± 2.3 


37.9 ± 1.8 




11.8 ± 2.4 


3.3 ± 0.7 


3.2 ± 0.7 


1.0 ± 0.2 


17.6 ± 3.5 


H 


Maffei2 


7.7 ± 1.6 


34.0 ± 1.4 




5.3 ± 1.1 


1.5 ± 0.3 


1.5 ± 0.4 


0.5 ± 0.1 


12.1 ± 2.4 


•2. 
a 


NGC1068 


14.0 ± 2.2 


38.5 ± 1.2 


11.5 ± 2.3 


16.0 ± 3.2 


10.4 ± 2.1 


3.9 ± 0.8 


0.9 ± 0.2 


21.3 ± 4.3 


NGC1097 


3.8 ± 0.9 


37.2 ± 1.7 




3.9 ± 0.8 




1.3 ± 0.3 




7.1 ± 1.4 




NGC1365 


7.0 ± 1.6 


35.4 ± 1.6 




5.5 ± 1.1 


2.7 ± 0.6 


2.3 ± 0.5 




11.0 ± 2.2 


(TJ 


NGC2146 


13.2 ± 3.8 


45.0 ± 2.6 


13.8 ± 3.1 


17.3 ± 3.5 


15.2 ± 3.1 


1.9 ± 0.4 


1.0 ± 0.2 


26.2 ± 5.2 


E. 


NGC253 


68.7 ± 12.0 


41.2 ± 1.4 


6.6 ± 1.6 


37.6 ± 7.5 


11.4 ± 2.5 


8.8 ± 1.8 


5.2 ± 1.1 


51.9 ± 10.4 




NGC3256 


7.1 ± 2.0 


42.0 ± 2.3 




12.8 ± 2.6 


6.0 ± 1.3 


1.3 ± 0.3 




13.7 ± 2.7 


S^ 


NGC3690 


8.0 ± 1.9 


51.3 ± 2.4 


7.0 ± 1.5 


8.3 ± 1.7 


8.5 ± 1.7 




0.5 ± 0.1 


8.6 ± 1.7 




NGC4038 


1.7 ± 0.6 


36.2 ± 2.7 




3.4 ± 0.7 


2.4 ± 0.5 






5.2 ± 1.0 


2- 


NGC4041 


1.7 ± 0.5 


35.1 ± 2.3 




1.9 ± 0.4 




0.5 ± 0.1 




3.5 ± 0.7 


^ 


NGC4414 


2.5 ± 0.4 


32.5 ± 1.0 




3.2 ± 0.6 




0.9 ± 0.2 




7.8 ± 1.6 


B 

(T) 


NGC4945 


46.2 ± 7.2 


36.1 ± 1.1 




19.3 ± 3.9 


7.5 ± 1.6 


6.9 ± 1.5 


3.4 ± 0.7 


35.2 ± 7.0 


& 


NGC520 


3.0 ± 0.8 


44.0 ± 2.4 




2.1 ± 0.4 


1.4 ± 0.3 


0.36 ± 0.08 


0.22 ± 0.05 


2.8 ± 0.6 


B 


NGC5430 


0.6 ± 0.4 


30.8 ± 3.9 




1.0 ± 0.2 


0.5 ± 0.2 


0.25 ± 0.06 




1.6 ± 0.3 




NGC5937 


0.8 ± 0.4 


32.4 ± 3.2 




1.5 ± 0.3 


0.8 ± 0.2 


0.33 ± 0.08 




4.0 ± 0.8 


<X> 


NGC6156 


1.3 ± 0.4 


32.4 ± 1.9 




1.8 ± 0.4 


1.1 ± 0.3 






4.1 ± 0.8 


^ 
^ 


NGC6240 


1.8 ± 1.4 


49.2 ± 7.7 




6.9 ± 1.4 




0.27 ± 0.07 


0.31 ± 0.07 


2.9 ± 0.6 


NGC6764 


0.3 ± 0.1 


29.8 ± 1.5 




1.1 ± 0.3 








1.1 ± 0.2 




NGC6810 


1.3 ± 0.5 


31.8 ± 2.4 




1.8 ± 0.4 




0.6 ± 0.2 




4.0 ± 0.8 


g 


NGC6824 


0.6 ± 0.2 


28.6 ± 2.5 




0.8 ± 0.2 




0.28 ± 0.06 




1.8 ± 0.4 


S' 

en 


NGC6946 


4.8 ± 1.1 


36.6 ± 1.6 




5.9 ± 1.2 


1.2 ± 0.3 


1.0 ± 0.3 


0.5 ± 0.1 


10.3 ± 2.1 




NGC7469 


1.6 ± 0.4 


40.8 ± 2.2 




2.5 ± 0.5 








2.3 ± 0.5 




NGC7552 


5.5 ± 2.0 


38.9 ± 2.8 




6.2 ± 1.2 


2.2 ± 0.5 






6.8 ± 1.4 




NGC7582 


3.3 ± 0.8 


37.1 ± 1.8 




3.0 ± 0.7 


2.7 ± 0.6 






4.2 ± 0.8 




NGC7714 


0.9 ± 0.5 


53.2 ± 5.8 




2.3 ± 0.5 


2.0 ± 0.5 






2.0 ± 0.4 




NGC891 


2.8 ± 0.9 


29.5 ± 1.9 




1.9 ± 0.4 




1.2 ± 0.3 




7.7 ± 1.5 





° Total far-infrared flux derived the dust temperature flt (see text). 

'' Dust temperature derived from the greybody flt (see text). 

" [OlIl]52/im line flux; "* [Ol]63/um line flux; " [OlIl]88/um line flux; ^ [Nll]122;um line flux; » [Ol]145/im line flux. '' [Cll]158/um line flux; 
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40 60 80 100 120 140 160 180 200 
Wavelength [|j.in] 

Fig. 1. The LWS spectrum and the fitted curve (Eq. (^) 
for M82. 

where Jl, T0.55, and B)^{T^) are the solid angle of the ob- 
ject, the optical depth at 0.55 /im, and the Planck func- 
tion of temperature Td, respectively. The parameters O 
and To. 55 cannot be determined independently and only 
the product of them is a meaningful parameter. We as- 
sumed that the dust emissivity is proportional to 1/A, 
which provides reasonably good fits for the present sam- 
ple of galaxies. A typical example of the fit is shown in 
Fig. ||. We then define the total infrared integrated flux 
from submicron dust grains as 



FIR^ 



F{X)dX. 



(2) 



Note that FIR does not include the excess emission 
shorter than 80 /im, which is attributed to the emission 
from very small dust grains. The uncertainties in the flux 
level of the SWl detector due to the memory effect and the 
spectral shape of the excess emission make it difficult to 
estimate the excess emission accurately. The excess emis- 
sion may be about 30% of FIR according to the study of 
the diffuse Galactic emission of Dwek et al. (1997). The 
present FIR is larger by 30% on average than the far- 
infrared flux for 42 — 122 /im estimated from the IRAS 60 
and 100 /zm flux densities (Helou et al. 1988). 

The observational results are summarized in Table g. 
In most of the sample galaxies, far-infrared forbidden 
lines, such as [Oi]63//m, [Oiii]88/xm, [Nil] 122 ^m, and 
[Cli]158/im, were detected. [Oi]145/im was detected only 
in 13 galaxies, while only 8 galaxies show detectable 
[O III] 52 fim emission. The indicated errors include the rel- 
ative errors estimated from the uncertainty in the base 
line fit and 20% flux uncertainty. The error in FIR is 
estimated from the uncertainty in the temperature deter- 
mination, which mostly comes from the goodness of the 
fit. 



[Cii]158/xm is one of the most important lines for the 
diagnosis of physical conditions of PDRs because of its 
large luminosity and low critical density for collisional ex- 
citation. However, the carbon atom has a lower ionization 
energy (11.26eV) than hydrogen (13.6eV), and carbon ions 
are expected to be present not only in the neutral region, 
such as PDRs, but also in the ionized regions. 

[Nil] 122 /im is a good tracer of diffuse low-density ion- 
ized gas, such as ELDWIM, because the ionization energy 
of nitrogen atom nearly equals to that of hydrogen atom 
and the critical density for collisional excitation by elec- 
trons is about 300 electrons cm^^. 

Oxygen atoms have an ionization energy of 13.6eV al- 
most the same as that of hydrogen. [Oi]63/im is one of 
the most luminous lines as well as [Cii]158/im and it be- 
comes a more efficient cooling line than [Cii]158/im in 
high-density gases. Together with the upper level transi- 
tion at 145 /im, it is an important probe for neutral gas. 
[O i]145/im was weak and detected only in a limited num- 
ber of galaxies, and we cannot examine the major frac- 
tion of the sample galaxies by using the [Oi]145/im line. 
[O III] 88 /im is a luminous line of dense ionized gas. It has 
an upper level transition at 52 /im and the line ratio of 
the 52 /im to 88 /im emission can be used to derive the 
electron density of the ionized gas (e.g., Moorwood et al. 
1980). Unfortunately, the spectra in the 52 /xm region do 



not have a sufficient signal to noise ratio to derive a re- 
liable [Oiii]52/xm line intensity for most of the present 
sample galaxies. 

In Table ra we list the electron density estimated from 
the ratio of the [O in] lines and the neutral hydrogen den- 
sity estimated from that of the [O i] lines for the galaxies in 
which [OIII]52 /im emission was detected. These are rough 
estimates and should be taken with caution because of the 
large errors in the obtained line ratios. For about a half of 
the galaxies with the detected [Oiii]52/im emission, the 
line ratio is near the low-density limit and only upper lim- 
its of the electron density are given. In the derivation of 
the neutral hydrogen density we assume that the gas tem- 
perature is lOOOK. Even with this temperature the line 
ratios are in the low-density limit for the galaxies listed 
in Table |^. For lower temperatures the upper limits will 
further be decreased. 

In Fig. |§i, the ratios of [Cii]158/xm and [Nil] 122 /xm 
fluxes to the far-infrared flux, [Cii]/F/i? and [Nu]/FIR, 
are plotted against the far-infrared color _R(60/100) while 
the ratio [Oi]63/im and [Oiii]88/im to FIR, [Oi]/ FIR 
and [O iii]/F/i?, are plotted against i?(60/100) in Fig. |b. 
Figure ||a indicates a trend that the ratio [C 11] /FIR de- 
creases with i?(60/100). The ratio [Nii]/F/i? also seems 
to decrease with i?(60/100). On the other hand, [O l]/FIR 
does not show a clear systematic trend with i?(60/100). 
There seems a weak trend that [O in] /FIR increases with 
the color, though the scatter is quite large. 

Similar trends in the ratios of the line intensities to the 
far-infrared intensity have bee n obtained fo r the normal 
galaxy sample (Malhotra et al. 1997 , 2001b ). The present 
sample includes not only normal galaxies but also star- 
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Table 3. Densities estimated from the line ratios 



object 



(cm 



^TTF" 



Cen A 


< 140 


< 56000 


Circinus 


170 ± 130 


< 10000 


M82 


200 ± 120 


< 1000 


NGC1068 


200 ± 120 


< 48000 


NGC2146 


125 ± 100 


< 45000 


NGC253 


<70 


c 


NGC3690 


< 180 


< 38000 



° The electron density derived from the [Olll]52/ini to 88 /im 

line ratio. 

'' The neutral hydrogen density derived from the [Ol]145/im 

to 63 fim line ratio for the gas temperature of lOOOK. 

'^ The line ratio of [O l] is too large and no reasonable density 

can be derived for NGC253 (see Fig. Rand next section). 



burst galaxies and AGNs. In Fig. ||, we plo t the data 
for normal galaxies by Malhotra et al. ( |2001b ) (filled cir- 
cles) together with the present sample (open squares) 
for comparison. In the figure, the total far-infrared flux 
F/i?[IRAS] is derived by the formula of Helou et al. ( |1988D 
based on the IRAS data to make direct comparison. Both 
show a similar trend to each other except that the data 
of normal galaxies seem to have a steeper trend with 
the color. AGNs (e.g., Cen A, Circinus, NGC1068, and 
NGC7582) are located in the same trend as normal and 
starburst galaxies, suggesting that the far-infrared emis- 
sion of AGNs is driven mainly by star-formation activities. 



3. Physical conditions of interstellar medium 

Together with the [Oi]63/im and 145 /xm lines 
[Cii]158^m can be used to derive the physical con- 
ditions of the line-emitting re gions based on PDR models 
(e.g., Tielens & Hollenbach 
HoUenbach & Tielens 



1997 



19851; Wolfire et 
Kaufman et al. 



al. 199C 



19991) , 



which the major model parameters are the incident FUV 
radiation field fiux Gq in units of t he so lar neighborhood 
value (1.6 x 10~^ W m~^, Habing 1968|) and the neutral 
hydrogen gas density n. However, the [C ii] line could 
also originate from ionized regions and the fraction of the 
contribution cannot be estimated a priori. 

To estimate the contribution to [Cii]158/im from 
PDRs, we take two approaches similar to Malhotra et 
al. (2001b). First we assume that all the emission of 
[Oi]63/im and far-infrared continuum (A > 80 /^m) comes 
from PDRs. Since the temperature of sub-micron sized 
dust grains is determined by the intensity of the incident 
radiation (e.g., Onaka 2000), Gq can be estimated from 
the dust temperature Td derived by Eq. (|l|). We used a 
semi-analytical equation of Go and Td given by HoUenbach 
et al. (1991) with Ay — 0.5 because a major fraction of 
[Cii]158/L(m and [Oi]63/im emis sions stem from the re- 
gion of Ay < 1 (Kaufman et al. 1999). 
rive Go = lO'^'^, while Kaufman et al. 



For M82 we de- 
estimated 



(119991) 




0.6 0.8 

R(60/100) 



1.2 



Fig. 2. a) Ratios of the [C ii] 158 ^m flux and [Nil] 122 ^m 
fiux to the total far-infrared flux FIR against the far- 
infrared color i?(60/100). Filled chcles indicate [C ii]/F/i? 
and open squares [Nii]/F/i?. b) Ratios of the [Oi]63/Ltm 
fiux and [Oiii]88/im flux to FIR against i?(60/100). 
Filled circles indicate [Oi]/ FIR and open squares 
[Oiii]/F/i?. 



tensities, suggesting that the present method provides a 
reasonable estimate of Gq. 

Then we compare the ratio of [O i]/FIR with the PDR 



model of Kaufman et al. (1999) with the derived Gq to es 



Go ~ 10 by taking account of several observed line in- 



timate n. Finally we estimate the intensity of [Cii]158/^m 
from PDRs with the derived Gq and n. The current esti- 
mate of the fiux uncertainty is 20% and there may be an 
uncertainty in the PDR model due to the assumed geom- 
etry. While FIR may be underestimated by a few tens % 
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Fig. 3. a) Ratio of the [Cii]158/im line flux to the total 
far-infrared flux FIR against the far-infrared i?(60/100). 
b) Ratio of the [Oi]63/im line flux to FIR against 
i?(60/100). Filled circles are taken from Malhotra et al. 
(2000) and open squares indicate the present sample. The 
value of i?(60/100) is based on IRAS data and FIR is 
derived by the formula of Helou et al. (1988). 



in the present analysis, it does not introduce a signiflcant 
error compared to other uncertainties. 

Figure ^ shows Gq and n against the color i?(60/100). 
The error in Go comes from the uncertainty in T^, while 
that in n is estimated from the uncertainty in [Oi] /i<"/i? 
to be about Alogn ~ 0.3 based on the FDR mod- 
els. Both parameters, Gq and n, increase clearly with 
i?(60/100) and n is found to be almost proportional to Go 
(logn = (0.8 ± 0.25) log Go + 0.75j:|^:4^, see also Fig. |l^. 
The trend mainly comes from the observed constancy 
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Fig. 4. a) FUV incident flux estimated by the dust tem- 
perature, Go, and the derived cloud gas density, n, against 
the far-infrared color i?(60/100). Filled circles indicate Gq 
and open squares n [cm~'^]. b) Ratio of the [Cll] intensity 
from PDRs to the observed total [C ii] intensity vs. the 
i?(60/100). 



of [Oi]/F/i? against the color. In the FDR model the 
photoelectric heating efficiency is roughly a function of 
Ga/n (Bakes & Tielens |1994D . For Go ex n, the heat- 
ing efflciency does not change and thus [Oi]/ FIR stays 
constant. The [Cii]158/im/F/i?, on the other hand, de- 
creases because the coUisional de-excitation becomes ef- 
ficient as n increases: the [CII] line has a lower critical 



density for excitation (~ 3 x 10 cm~^) than [Oi]63/^m 
(~ 1 X 10^ cm""^). In this situation, we expect an increase 
in the gas temperature. Because the gas temperature is al- 
ready sufficiently high in the parameter range in question, 
the increase in the temperature does not affect the line 
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intensities appreciably. Hence the decrease in [Cii]/F/i? 
with color can be attributed to the increase in the coUi- 
sional de-excitation of the [C ii] transition for the present 
sample of the galaxies. The [Cii]158/im emission from 
PDRs is estimated to be '^ 50 ± 30% of the total observed 
[Cil]158/zm emission in the present analysis and the frac- 
tion of the PDR component is indicated to decrease with 
i?(60/100) (Fig. |b). Therefore, the decrease in the total 
[C ii] /FIR is attributed to the decrease in the PDR com- 
ponent due to the thermalization of the [C ii] transition. 

In order to examine the reliability of the present data 
reduction and analysis we can compare the present results 
with those of previous works for some individual galaxies. 
The comparison is summ arized in Table IJ. Colbert et al. 
(1999) and Unger et al. (|200C| ) analyzed the same LWS 
data and derived Go and n by using the PDR models 
of Kaufman et al. (1999) for M82 and Cen A, respec- 
tively. The present analysis provides fairly good agree- 
ment with their results. The line intensities derived in 
the pres ent s tudy also agree with those by Braine & 
Hu ghes (1999 ) within the measurement errors. Carral et 
al. ( 1994 ) reported the results of observations of NGC253 
and NGC3256 by the KAO and obtained Go and n based 
on the PDR models of Wolfire et al. (|1990D. While the hue 



fluxes are in agreement with the present results within the 
errors except for [Oiii]88/xm of NGC253, the present re- 
sults indicate systematically low densities. The difference 
can be attributed to the relatively high gas temperature 
with the same Go and n in the models of Kaufman et al. 
because their models include the additional gas heating 
due to polycyclic aromatic hydrocarbons (PAHs). 

The [Cii] emission other than the PDR origin may 
be ascribed to the extended low density warm ionized 
medium (ELDWIM). A large fraction of [Nil] 122 ^m hue 
emission is tho ught to stem most ly from the ELD WIM 
(Wright et al. |l99ll He iles |l994t Bennett et al. |l994 ; 
Petuchowski et al. 1994 ). Figure |^ plots the ratio of the 
observed [Nil] 122 /im intensity to the non-PDR compo- 
nent of [Cii] 158 /im. The ratio shows a large scatter with 
a weak trend that the ratio decreases with i?(60/100). 
Because of the uncertainties the trend may be spurious 
(see below). 

All the data points are located within the range 0.1 — 
0.7 in Fig.|5[ The line ratio of [Nii]/[Cii] in the ionized gas 
depends on the electron density, but is insensitive to the 
temperature of the ionized gas. To estimate the line ratio 
expected from ionized gas we assume To = 7000K and the 
abundance of [C+]/[H+] = 4 x l{)-^5c+ and [N+]/[H+] 
= 1 X 10^^(5n+ with the depletion factors as 5q+ — 0.65 
and (5n+ = 1 in the following discussion (Heiles 0994). 
We adopted the collision coefhcients for N+ from Stafford 
et al. ( |1994| ) and those for C+ from Heiles (|l99§). Recent 
-/KT observations indicate that the interstellar abundance 
of carbon and nitrogen in the gas phase is fairly constant 
on various l ines o f sight in our Galaxy as ^c = 0.35 ±0.0 5 
(Sofia et al. |1997|) and 5^ = 0.75±0.05 (Meyer et al. |l997D . 
Based on these values the relative abundance of N to C 
will increase by 30%. The following discussion thus has 
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Fig. 5. Ratio of [Nil] 122 /im intensity to the non-PDR 
component of [Cii]158/im against i?(60/100). 



an uncertainty of this level associated with the assumed 
abundance. 

The lower boundary of the observed ratio 0.1 is then 
found to correspond to the low-density limit of the ratio 
in the ionized gas. The upper boundary 0.7 is obtained 
for a gas with Uc — 120 cm^'^. Petuchowski et al. (1994) 
reported a large [Nii]122/im to 205 /im line ratio in the 
central 850pc of M82 compared to the average ratio of the 
Milky Way (Wright et al. 1991), indicating that a fair frac- 
tion of the [N ii] line emission comes from the ionized gas 
of Tie = 150- 180cm-3 in M82. The ratio of [Nil] 122 /im 
to [C ii]158/im of non-PDR origin for M82 is about 0.3 in 
the present analysis, suggesting that there may be a signif- 
icant contribution to the [CII] emission from the low den- 
sity diffuse ionized gas in the outer part (> 850pc) of the 
galaxy. The observed intensity is compatible with the in- 
terpretation that the non-PDR component of [Cii]158/im 
comes from the ionized gas that emits [Nil] 122 /tm. 

The ionized gas also emits a radio continuum. The in- 
tensity of free-free transition is written for T^ — 7000K 

by 



/(4.85GHz)[mJy sr" 



3.12 X IQ-^^nll 



(3) 



where ric is the electron density in cm~^ and I is the path 
length in cm (Spitzer 1978). The intensity of the [Cii] line 
from the ionized gas is given by 



/r,^ jji (ELDWIM) = {l/^-K)L{T^)n(C+)nJ, 



(4) 



where L{Tc) is the co oling function of [Cii]158/im (e.g., 
Hayes & Nussbaumer 1984). For a given [Cii]158/im in- 
tensity, the corresponding radio continuum intensity in- 
creases with ric as can be estimated through Eqs. (y) and 
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Table 4. Comparison with previous studies 



object ref." [Ol]63^m [Oni]88^m [Nll]122^m [Ol]145/im 
(IQ-'^^Wm-^) (IQ-'^Wm-') (IQ-'^Wm"') (lO'^^Wm-^; 



"C II] 158 Mm Go 



(10"'^Wm 



(cm 



M82 


P 


169 ± 34 


91 ± 18 


21 ± 4 


15 ± 3 


128 ± 26 




1 


176 ± 5 


86 ±4 


17 ± 3 


12 ± 1 


134 ± 1 


Cen A 


P 


19 ± 4 


6.3 ± 1.4 




1.0 ± 0.2 


29 ± 6 




2 


19.6 


7.0 


1.5 


1.1 


29.1 


NGC4414 


P 


3.2 ± 0.6 




0.9 ± 0.2 




7.8 ± 1.6 




3 


3.3 ± 1 


1 ± 0.5 


1.3 ± 0.4 




10.6 ± 2 


NGC253 


P 


38 ± 7.5 


11.4 ± 2.5 






52 ± 10 




4 


45 ± 6 


6± 1 






48 ± 2 


NGC3256 


P 


12.8 ± 2.6 


6.0 ± 1.3 






13.7 ± 2.7 




4 


14.3 ± 2.6 


4.8 ± 1.2 






11.7 ± 2.5 



10^ 
10^ 



10^ 



lO'^ 

lO-*' 
10^- 
10^ 



lO'' 
10^ 
10^ 



10'^ 



10"' 
10* 
lO^' 
10^' 



references: p present work; 1 Colbert et al. (1999); 2 Unger et al. (200C); 3 Braine & Hughes (1999); 4 Carral et al. (1994) 
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Fig. 6. The observed 4.85GHz radio continuum flux 
density plotted against the 4.85GHz radio continuum 
flux density predicted from the non-PDR component of 
[Cii]158/im emission (see text). The observed radio flux 
densities are taken from Gregory et al. (1991, 1994), 
Becker et al. (1991), and Griffith et al. (1994). 



(0). As a conservative upper limit we assume that the non- 
PDR component of [C ii] 158 /im emission comes from the 
ionized gas of density 200 cm~^. Figure |6| shows the com- 
parison between the predicted and observed flux densities. 
We take the 4.85GHz radio continuum data from Becker 
et al. ( |1991D , Gregory & C ondon ( |1991| ), Gregory et al. 
( |1994D , and Griflith et al. ( |1994| ). The beam size of the 
radio observations was about 4' (FWHM). Observations 
of 1.49 GHz and 1.425 GHz indicate that most of the ra- 
dio continuum emission comes from the area within about 
1^ of the c enter of the infrared emission (Condon et al. 



199C, 1996 ). We assume that the observed 4.85 GHz emis- 
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Fig. 7. Ratio of the intensity of the [Oi]145/im line to 
the [Oi]63Mm line against i?(60/100). 



The predicted value shows a positive correlation with the 
observed flux density except for Ccn A. Cen A is a very 
strong radio source and most of the radio emission from 
Cen A is nonthermal (e.g. Sreekumar et al. 199£). The 



sion also comes from the region within the LWS beam. 



observed radio intensities are larger than the upper lim- 
its predicted from [Cii]158/xm emission of the non-PDR 
origin for most of the sample galaxies. Thus the non-PDR 
component of [C ii] 158 /xm emission is compatible with the 
observed radio continuum intensity when it arises mostly 
from the low-density ionized gas for the present sample of 
galaxies. 

Figure M shows the ratio of [Oi]145/im to [Oi]63/im 
against i?(60/100). The observed ratio is in the range 
0.05 — 0.1 for most galaxies, which is in agreement with 
the prediction of the PDR model with the range of Go 
and n estimated in the present analysis. The ratio does 
not exceed 0.1 within a re asona ble range of the density 
and temperature (Watson 1984 ). However, a few galax- 
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ies clearly show higher ratios than the model. They are 
NGC4945 and NGC253, both of which have large incli- 
nation angles (for NGC4945 and NGC253, i ~ 85° and 
80°, respectively). NGC520, which shows the third largest 
ratio, has i = 60°. The absorption of [Oi]63/ini in the 
interstellar medium may affect the ratio in these galax- 
ies. For a majority of the galaxies, however, this effect is 
probably not significant. 

In the analysis described above it is difficult to prop- 
erly evaluate the uncertainties in the comparison with the 
model. In order to examine how robust the derived conclu- 
sions are, we take another approach to estimate the PDR 
contribution to the [Cii]f58/im line emission. Figure |5| in- 
dicates that there seems no strong trend in the ratio of 
[Nil] 122 /^m to the non-PDR origin of [Cii]158^m. We 
thus simply assume that the contribution from the ion- 
ized gas to [Cii]158/^m emission is proportional to the 
[Nil] 122 /im intensity. We take a mean value of Fig. ^ as 



[Cii]158Afm(ELDWIM) = 3.5 x [Nil] 122 /im. 



(5) 



This relation corresponds to a gas of ric ~ 35cm ^ for 
Te = 7000 K. We estimate the intensity of [Cii]158/im 
from the ionized gas based on the [Nil] 122 /im intensity 
by using Eq. (bh and attribute the rest to that coming 
from PDRs. We then estimate Gq and n from [Oi]/[Cii] 
and ([Cii]-f [Oi])/i^/i?. Figure ||a shows Gq and n esti- 
mated in this method. The values of Gq derived in this 
method are in agreement with those obtained in the first 
approach within the estimated errors. Thus the param- 
eters derived in the second method are still compatible 
with the continuum spectrum of the LWS spectra. We 
have obtained the same trend as in the first approach: 
both n and Go increase with the color and n is roughly 
proportional to Go (see also Fig. ^0|). We conclude that 
the linear increase of n with Gq is a rather secure result 
for the present sample galaxies. The weak trend seen in 
Fig. |5| is not necessarily real. Figure gb shows the frac- 
tion of the PDR component of [C ii] emission derived in 
this analysis, suggesting that the PDR component does 
not vary with _R(60/100) in contrast to Fig. |jb. Thus in 
this analysis the decrease in [C ii] /FIR can be interpreted 
in terms mainly of the decrease in the ionized component 
relative to FIR as indicated in the decrease in [N Ii]/F Ji? 
(Fig. la), though the decrease in [Cii](PDR)/F/i? due to 
the thermalization also contributes partly. 

The ratio of [Cli]158/im to ^^CO (J=l-0) line inten- 
sity is another measure for the diagnosis of PDRs. Because 
C^ is converted to CO as the gas shielding to prevent CO 
from dissociation becomes efficient, the thickness of C^ 



layer is a function of Gq/u (Mochizuki & Onaka 2001). 
The intensity ratio increases with Go and decreases with 
n and it is a function of Go/n in the range of Go and 
n in question (Pierini et al. 1999; Kaufman et al. 1993; 



Mochizuki & Nakagawa |2000|) 
impor tant role in the C 
et al 



Metallicity also plays an 



to CO conversion (Mochizuki 
1994). Hence, unless there is an appreciable vari- 
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Fig. 8. The results of the analysis based on the as- 
sumption that the contribution from the ionized gas to 
[Cii]158/im is proportional to [Nil] 122 /im as given by 
Eq. d) (see text), a) FUV incident flux. Go (filled cir- 
cle) and the gas density, n (open square) derived from 
[Oi]/[Cii] and ([Cii]-t-[Oi])/F/i? of the PDR model by 
Kaufman et al. (1999) against i?(60/100). b) The ratio of 
the estimated [C ii] intensity from PDRs to the observed 
total [Cii] intensity against i?(60/100). 



plots the ratio against the color. The CO data are taken 
fro m You ng et al. ( |1995D , Elfhag et al. ( |1996D , A aho e t 
al. ( p^95| ), Curran et a l. ( pOOll ), Stacey et al. ( p^9l| ), 
and Eckart et al. (1990). The CO observations had the 
beam size of 45 



ation in the metallicity in the sample galaxies, the ra- 
tio of [C ii]/CO approximately varies with Gq/ti. Figure ^ 



56". We assume that the CO emission 
comes mostly from the central part of the galaxies and 
did not make any corrections for the beam size. Except 
for NGC6824, the ratio stays almost constant, supporting 
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Fig. 9. Ratio of the [Cii]158/im intensity to the CO 
(J=l-0) intensity against i?(60/100). The CO emission has 
not been detected for NGC6824 and a lower hmit for the 
ratio is plotted. 

n oc Gq. The constancy of Gq/tl is also confirmed by the 
[Cii]/CO ratio. The CO (J= 1-0) emission in NGC6824 
was not detected (Young et al. [1995 ). With the upper 
limit the ratio of [Cii]/CO is estimated to be larger than 
10^, which is in a simila r rang e to those found in quiet spi- 
rals (Smith & Madden 1997 ) . The values of the observed 



ratio of the other galaxies 8 x 10 — 3 x 10 are in agree- 
ment with those observed in Galactic PDRs (Stacey et al. 
1991). It is slightly smaller than the values predicted in 



the PDR model. 



4. Discussion 



Malhotra et al. ( |1997| , |2001b ) have examined six possi- 
ble interpretations for the observed decrease in [Cii]/F/i? 
with the FIR color _R(60/100) and favored the interpre- 
tation that the charge up of dust grains reduces the ef- 
ficiency of photoelectric heating. When dust grains are 
charged positively under strong ultraviolet radiation, the 
Coulomb potential prevents electrons from escaping and 
thus the photoelectric heating efficiency decreases. The 
electric charge of dust grains is determined by the bal- 
ance between ionization and recombination and Go/n is a 
good measure to indicate the electrical potential of dust 
grains and thus the heating efficiency in the parameter 
range obtained in the present study, Gq ~ 10^ — 10'' and 
n ~ 10^ - lO^cm-^ (Bakes 
1999| ). 



Tielens 1994; Kaufman et 



To examine this hypothesis, we plot Go/n against 
i?(60/100) in Fig. |o|. There is no clear trend seen in 
the figure, suggesting that at least for the present sam- 
ple of galaxies Gq/ti does not increase with the color. The 
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Fig. 10. The ratio of Go to n against the far-infrared color 
i?(60/100). Filled circles indicate the results derived by 
the comparison of [O i] /FIR with the PDR model, while 
open squares depict those derived by assuming that the 
[Cii] from the ionized gas is proportional to [Nil] 122 /im 
emission (Eq. (||)). 



present sample contains starbur st gala xies compared to 
the sample of Malhotra et al. ( ^OOlb ). The charge up 
of dust grains is supposed to occur more easily in star- 
bursts than in normal galaxies. However, the decrease 
of [Cn]/ FIR is more steep in Malhotra et al. sample 
than in the present sample (Fig. ^ja). The interpretation 
of charge up of dust grains seems to not play a signifi- 
cant role in the trend of [Cu]/FIR. Instead we suggest 
that the decrease could be attributed to either the in- 
crease in the de-excitation of the [Cii]158/im and/or the 
decrease in the ionized gas component. The former pos- 
sibility has been pointed out also by Genzel & Cesarsky 
( 2000| ) . The latter possibility can be examined by observa- 
tions of [Nii]205^m which together with the [Nil] 122 /^m 
intensity allow a better estimate of [Cii]158/.tm emission 
from the ionized gas. The decrease in the ionized gas com- 
ponent may be a consequence of the relative increase in 
the far-infrared intensity related to the activities. Helou 
et al. ( ^001 ) have also suggested that the decrease of the 
ratio could result from the decrease in the importance of 
PAHs or small grains in gas heating based on the observed 
correlation between [Cii]158/im and the 5 — 10 /im inten- 
sities. 

Based on the observations of Ha, Hi, CO, and far- 



infrared continuum, Kcnnicutt (1998) has investigated the 
global Schmidt law in galaxies and found that the disk- 
averaged star formation surface density Ssfr is propor- 
tional to S^g with N — 1.4 ± 0.15, where Sgas is the 
surface density of the gas. The present results indicate 
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Go oc n^ with M = 1.25+°:3 (see Sect. |). If Go simply 
indicates the global star-formation rate (SFR) in galaxies, 
this relation suggests a similar dependence of SFR on the 
volume gas density of star-forming regions and little vari- 
ation in the scale height of star-forming gas distribution 
in galaxies. 

The present galaxy sample contains AGNs, star burst, 
and normal galaxies, and the present results are quite sim- 
ilar to those for the normal galaxy sample (Malhotra et al. 
2001b ). Thus they are thought to indicate general char- 
acteristics of far-infrared properties for a wide range of 
galaxies, though irregular, early- type, and quiescent spi- 
ral galaxies may s how s lightly differe nt ch aracteristics (see 
Smith & Madden |l997|; L eech et al. |1999| ; Malhotra et al. 
2001a ; Hunter et al. 2001 ). We do not see any clear differ- 
ence in the trends for AGNs in the present sample, sug- 
gesting that the far-infrared properties in AGNs are also 
driven by star-forming activities. A similar conclusion has 
been drawn for Cen A by Unger et al. ( |2000| ). ISOPHT 
observ ations of CfA Seyfert sample by Perez Garcia et al. 
( 1998 ) also indicate that FIR SED of Seyferts can be inter- 
preted in terms of the thermal emission from star-forming 
regions. 

We summarize the general trends in the far-infrared 
properties: 

[Gii\/FIR = (l.ltoi) X 10-3i?(60/100)-(°-^^-o'2), (6) 

[Oi]/F/i?= (1.25 ±0.75) X 10-^ (7) 

[Nii]/F/i?= (l-2to:5) X 10"''i?(60/100)~(^-^-o*), (8) 



and 



[Om]/ FIR = (l.i: 



0.8/ 



10-3i?(60/100)(i'^-°?s). (9) 



The last relation of [Oiii]/F/i? has a large scatter and 
should be taken with caution. 

The present analysis also suggests a general relation 
between the typical density in the galaxy and the far- 
infrared color as 



(3.0t? 



+2.0 



5; X 10-^i?(60/100) 



(2.4±0.6) 



(10) 



These relations suggest a possibility that the physical 
properties of distant galaxies can be estimated solely from 
the far-infrared color. They can also be used in the com- 
parison of the physical properties of distant galaxies to 
examine differences from nearby galaxies in the activities, 
if any, in observations by future space missions, such as 
SIRTF, ASTRO-F, and Herschel Space Observatory. 

5. Summary 

We have investigated the LWS spectra of 43 — 197 /im 
for 34 nearby galaxies. In addition to the detected emis- 
sion line intensities, we estimated the far-infrared inten- 
sity from submicron grains from the continuum emission 
of > 80 /im and derived the average dust temperature. We 
obtained the following results. 



We have found that the ratios of [Cil]158//m and 
[Nil] 122 /im flux to the total far-infrared flux FIR 
decrease as the far-infrared color i?(60/100) becomes 
bluer, but the ratio of [Oi]63/im to FIR does not 
show a systematic trend with the color. The ratio 
of [Oiii]88/im to FIR shows a large scatter with a 
weak trend of increase with the color. These are simi- 
lar to those obtained for the normal galaxy sample by 



Malhotra et al. ( pOOlbj ). 

2. There is no clear difference seen between AGNs and 
starburst galaxies in the present sample in these 
trends, suggesting that even in AGNs the far-infrared 
properties are dominated by star-formation activities. 

3. Based on the c ompa rison with the PDR model by 
Kaufman et al. ( 199£ ), we found that the typical neu- 
tral gas density in galaxies increases linearly with the 
radiation field intensity. 

4. About a half of the observed [Cii]158/im emission is 
estimated to come from PDRs. The rest can be at- 
tributed to that coming from the diffuse ionized gas 
that emits the [Nil] 122 /im line. 

5. The observed decrease in [Cii]158/im/i^/i? with the 
FIR color can be interpreted in terms of either the in- 
crease in the collisional de-excitation in the [C 11] tran- 
sition at high density or the decrease in the ionized 
component. Decrease in the photoelectric yield due to 
the charge up of dust grains does not seem to play a 
primary role in the observed trend. 

6. We summarize the relations among the far-infrared 
properties for the present sample galaxies. These are 
thought to indicate the general characteristics for a 
wide range of galaxies, including starburst, AGNs, and 
normal galaxies and can be used as a measure in the 
investigation of activities in distant galaxies in future 
observations. 
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